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ABSTRACT 

Context. To form metre-sized pre-planetesimals in protoplanetary discs, growing grains have to decouple from the gas before they are 
accreted onto the central star during their phase of fast radial migration and thus overcome the so-called "radial-drift barrier" (often 
inaccurately referred to as the "metre-size barrier"). 

Aims. To predict the outcome of the radial motion of dust grains in protoplanetary discs whose surface density and temperature follow 
power-law profiles, with exponent p and q respectively. We investigate both the Epstein and the Stokes drag regimes which govern 
the motion of the dust. 

Methods. We analytically integrate the equations of motion obtained from perturbation analysis. We compare these results with those 
from direct numerical integration of the equations of motion. Then, using data from observed discs, we predict the fate of dust grains 
in real discs. 

Results. When a dust grain reaches the inner regions of the disc, the acceleration due to the increase of the pressure gradient is 
counterbalanced by the increase of the gas drag. We find that most grains in the Epstein (resp. the Stokes) regime survive their radial 
migration if -p + q + j < (resp. if q < |). The majority of observed discs satisfies both -p + q + ^ < and g < | : a large fraction 
of both their small and large grains remain in the disc, for them the radial drift barrier does not exist. 

Key words, planetary systems: protoplanetary discs — methods: analytical 



■ 1. Introduction This migration motion depends strongly on the grain size, 
QQ I which sets the magnitude of the drag, a s well as the nature of the 
O ' Much of the information about the gas structure of protoplane- '^'"^8 '"^8™^- Specifically, as shown by^ and^iSHS^ grains 
^. tary discs is inferred from the emission by the dust component ^ '^"^^'^^1 ^'^^ P^^^ through the disc m a fraction of the disc 
^ ■ and an assumed dust-to-gas ratio. Interpr etations of recent ob- J^f^^^f,- T^is catastrophic outcome is ca led the radial-drift 
^ • servations i n the (sub)minim etre domain (lAndrews & WiUiamsl b^™^"" °f Pl^"^^ formation. More precisely, we will adopt the 
I2OOI l2007t lLommen etal.ll2007h show that observed discs typ- subsequent definition for the 'radial-drift barrier in this study: 
ically have masses between 10-^ and lO"' and a spatial 'the ability of grains of be accreted onto the central star/deplete 
extent of a few hundred AU. Their radial surface density and ^0"^ '^'^'^ lif^t™^ ■ Historically, this process 
• , temperature profiles are approximated by power laws (2 cc r", ^''^^ ^^^'^i^'l ^ Minimum Mass Solar Nebula (MMSN,^ 
X • r ^ whose respecliJe exponents /and q have positive s££^eidenschilling 1977b; Havashi 1981; Desch 2007; CridJ 
b : values typically of order unity. l2QQ9|), in which the critical size corresponds to metre-sized bod- 
5t 1 ies and thus was called the "metre-size barrier". However, plan- 
Seminal studies describe the dust motion in protoplan- ets are frequently observed (besides the 8 planets in our solar 
etary discs, wh i ch dep ends strongly on th e gas structure system more than 700 extra-solar planets have been discovered 
IWeidenschillingI (Il977ai hereafter W77) and iNakagawa etall to datfl): some solid material must therefore have overcome this 
19861 hereafter NSH86') demonstrated that dust grains from mi- barrier and stayed in the disc to form larger bodies. Moreover, 



cron sizes to pre-planetesimals (a few metres in size) experi- if the small grains of every disc were submitted to the radial- 
ence a radial motion through protoplanetary discs. This motion drift barrier, we would barely detect them since their emission 
is called radial drift or migration. Due to its pressure gradient, via optical/IR scattering and IR thermal radiation is due to small 
the gas orbits the central star at a sub-Keplerian velocity. Grains grains. As discs are frequently observed, the grains from which 
therefore have a diff^erential velocity with respect to the gas. The the emission is detected cannot be strongly depleted for a sub- 
ensuing drag transfers linear and angular momentum from the stantial fraction of discs, 
dust to the gas. Thus, dust particles can not sustain the Keplerian 

motion they would have without the presence of gas and as a 

result migrate toward the central star. ' http://exoplanet.eu 
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From a theoretical point of view, such a discrepancy be- 
tween the observations and the theoretical predictions imply that 
the seminal theory has to be extended (some physical element 
is lacking) or that it has not be en fully exploi t ed. Th is second 
option has been investigated by lYoudin & Shul (l2002l hereafter 
YS02). Th ey high light the fact that, contrary to the primary hy- 
pothesis of |W77L observed dusty discs are drastically different 
from the MMSN prototype. As the radial surface density and 
temperature profiles fix both the radial pressure gradient and the 
magnitude of the gas drag, different values for the power-law 
exponents p and q affect the optimal grain size of migration and 
thus induce different radial motions for the dust through the disc. 
Specifically, YS02 showed that for steep surface density profiles 
and smooth temperature profiles, the grains radial velocity de- 
creases when the grains reach the inner discs regions. Grains in 
such discs therefore experi ence a "pile-up". However, while im- 
portant, the work of lYS02l does not provide a precise conclusion 
on the outcome of the grains nor any quantitative criterion for the 
"pile-up" process to be effic ient enough to avoid the radial-drift 
barrier Furthermore, |YS02| restricts their study to the special 
case of a gas phase with a low density (e.g. the grain size smaller 
than the gas mean free path, called the Epstein regime). This hy- 
pothesis is not valid anymore when considering the radial drift of 
pre-planetesimals, whose grain sizes are larger than the gas mean 
free path and are submitted to the Stokes drag regime. Although 
the radial drift of pre-planetesimals has already been studied in 
different s ituations with numerical or s e mi-analytical method s 
— see e.g.lHaghighipour & BossI (l2003l) : lBirnstiel et alJ (l2009h : 
lYoudrnl (l2011h — its rigorous theory for the standard case of a 
simple disc has not yet been derived. 

Within this context, we see that (i) the seminal theory de- 
scribing the radial motion of dust grains has been developed 
within the limits of the Epstein regime but does not treat the 
Stokes regime ; (ii) here exists no clear theoretical criterion to 
predict the impact of the "pile-up effect" on the outcome of 
the dust radial motion ; (iii) there exists no criterion to predict 
whether a given disc will be submitted to the "radial-drift bar- 
rier" phenomenon. To answer these three points, we re-visit in 
this study t he wor k of IW77l and NSH86 and extend the devel- 
opments of IYSO2I for both the Epstein and the Stokes regime. 
Performing rigorous perturbative expansions, we find two the- 
oretical criteria (one for each regime) which predict when the 
"pile-up" effect is sufficient for the grains not to be accreted onto 
the central star. We then test when these theoretical criteria can 
be applied in real discs. 

Additionally , our work i s motivat ed by the recent observa- 
tional results of iRicci et al ] (I2010albh . From their observations 
they claim that "a mechanism halting or slowing down the in- 
ward radial drift of solid particles is required to explain the data". 
In this work we aim to show that contrary to what is usually in- 
voked, local pressure maxima due to turbulent vortices or spi- 
ral density waves may h elp but are not nece ssarily required to 
explain the observations. iRicci et a n (I2010ah also mention that 
"the observed flux of the fainter discs are instead typically over- 
predicted even by more than an order of magnitude". Here, we 
also aim to provide a quantitative criterion to determine which 
discs are faint and which one are not. Thus, revisiting the sem- 
inal theory of the radial drift is timely, all the more so than an 
important quantity of new data is soon to be provided by ALMA, 
the Atacama Large Millimeter/submillimeter Array. 

In this paper, we first recall some general properties of grain 
motion in protoplanetary discs for both the Epstein and Stokes 
regime in Sect. |2] We then focus on the radial motion of non- 
growing grains in the Epstein regime. We expand the radial mo- 



tion equations assuming a weak pressure gradient in Sect. [3] and 
detail the two different modes of migration which grains may 
experience in Sects. l3.1l and [T2l This allows us to derive an an- 
alytic criterion which determines the asymptotic dust behaviour 
in the Epstein regime in Sect. 13. 31 We transpose these derivations 
for the Stokes regime at low Reynolds number in Sect.|4]and ob- 
tain a similar criterion for this regime. We also discuss the grains 
outcome for large Reynolds numbers. In Sect.|5] we discuss the 
relevance of these criteria and study their implications for ob- 
served protoplanetary discs and planet formation in Sect.|6] Our 
conclusions are presented in Sect. |7] 



2. Dynamics of dust grains 

To reduce the parameter space for this study, we assume the fol- 
lowing: 



1. The disc is a thin, non-magnetic, non-self-graviting, invis- 
cid perfect gas disc which is vertically isothermal. Its radial 
surface density and temperature are described by power-law 
profiles. Notations are described in AppendixlAl The flow is 
laminar and in stationary equilibrium. Consequently, the gas 
velocity and density are described by well-known relations, 
which we present in AppendixlB] 

2. The grains are compact homogeneous spheres of fixed ra- 
dius. The collisions between grains and the collective effects 
due to large dust concentrations are neglected. When the 
grains are small compared to the mean free path of the gas 
(/Ig > 4i/9, where s is the grain size), their interactions with 
the gas are treated by the Epstein drag force for diluted me- 
dia (E pstein 1924; Baines et al. 1965; Stepinski & Valageai 
Il996l) . This drag is caused by the transfer of momentum by 
individual collisions with gas molecules at the grains sur- 
face. Assuming specular reflections on the grain and when 
the differential velocity with the gas is negligible compared 
to the gas sound speed, the now common expression of the 
drag force is 



Fd Av 



PgCs 



(1) 



where md is the dust grain's mass, the stopping time, pg 
the gas density, Cj the local gas sound speed, pd the intrin- 
sic dust density, and Av = v - Vg the differential velocity 
between dust and the mean gas motion. In classical T Tauri 
star (CTTS) protoplanetary discs, drag forces for particles 
smaller than ~ 10 m are well described by the Epstein regime 
( iGaraud et al.ll2004l see also Sect. 16.1b . Small grains which 
produce the emission of observed protoplanetary discs sat- 
isfy this criterion. 

The interactions between large dust particles (As < 4s /9) and 
the gas are treated by the S tokes drag force dWhippldl 19721 : 
IStepinski & Valageaslll996l) . In this case, the gas mean free 
path is small and the dust particle is locally surrounded by 
a viscous fluid. Depending on the local Reynolds number 
of the flow around the grains Rg - where v is the 

microscopic kinematic viscosity of the gas, the drag force 
takes the following expression: 



Fd = --CDns-pg\A\\ Av, 



(2) 
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where the drag coefficient Co is given by 
24/?:' for/?, <1 



Cr 



2ARf -^ for 1 < /?g < 800 



(3) 



0.44 for 800 < . 
\f Rg< I, the drag force remains Hnear in Av. 

In this work, the physical relations are written in cylindrical 
coordinates (r, 9, z). The related unit vector system is given by 
(e,., Ce, Cj). As the system is invariant by rotation around the ver- 
tical axis C;, the physical quantities depend only on r and z- The 
physical quantities of the gas, designated by subscript g, are first 
determined in a general way. Then, the limit z = is taken to 
study the restricted radial motion. 

Dust dynamics depends on both the magnitude of the drag 
(driven by the differential velocity) and on its relative contri- 
bution with respect to the gravity of the cen tral star S emina l 
studies of dust dynamics were conducted by IWhippld (Il972h . 



ers dYSOl Takeuchi & Lin 20021 Ha 


shishiDOur & Bossll2003l 


iGaraud et al. l2004l 1 Youdin & Chiand 


20041). Here we recall the 



major points of those studies. We consider two forces acting on 
the grain: the gravity of the central star and gas drag. (We as- 
sume that the momentum transferred by drag from a single grain 
on the gas phase is negligible.) Thus 



'"d— = -Fd +'«dg, 
df 



(4) 



where Fd is the drag force. As shown by Eqs. ([T)-©, a general 



expression of the ratio ^ is of the form 



Fd 



C(r,z) 



V - Vgl-* (v - Vg) , 



(5) 



where the quantities C, y and A are defined for both the Epstein 
and the Stokes regime in AppendixICl 

3. Radial motion in the Epstein regime: 

perturbation analysis at small pressure gradients 

Considering the Epstein (small grains) regime, Eq. (|4| reduces 
to: 

OTd— = V - Vg -H mdg. (6) 

df fs ^ ^' 

Writing Eq. ^ in (r, 0, z) coordinates leads to 
{\'r - Vgr) gMr 



dVr 

"d7 



r 



"d7 



df 



(z2 + r^f- 

VrVg (vfl-Vgfl) 



+ 



(7) 



fs 

SMz 



(z2 + r2) 



3/2 



To highlight the important parameters involved in the grains dy- 
namics, we introduce dimensionless quantities (see Appendix 
let . It is crucial to note that the ratio ^ of the two timescales 
related to the physical processes acting on the grain is given by 

- = — /J^ei^ = SqRPs^ . (8) 



With Eq. ([T]) and noting Qk the Keplerian angular velocity, this 
ratio can be written as 



fs 



(9) 



^opt 



where s, 



opt 



This timescale ratio therefore corresponds to a 



dimensionless size S - S^R^e^K^-'i for the grain. If 5 «: 1 (resp. 
5 » 1), the effects of drag will occur much faster (resp. slower) 
than gravitational effects. If S - 1, both gravity and drag will 



act on the same timescale. Interestingly, Sopi 
midplane as r^P, as does surface density. 

Then using Eqs. ©, dRlb and dCTT ) 
(er,efl,ej): 



dfv Vg Vr /„j-^\ z- R 



varies in the disc 



dr R So 

dVfl VgVr 
1 + 

dr R 



(/?2 + 0oZ2) 



we obtain for 



= 



3/2 



(10) 



d^Z 1 dZ /„_,3N Z 

5odr {R^~ + ^oZ^f' 




0. 



These equations depend on five control parameters: the initial 
dimensionless grain size, ^o, the radial surface density and tem- 
perature exponents, p and q, the square of the disc aspect ra- 
tio, (po - (//(j/ro)^, and the subkeplerian parameter, 770, given by 
Eq. iC.2i . The equations can be simplified in some cases, e.g. if 
the vertical motion is considered to occur faster than the radial 
motion, R ^ 1 and ^ simplifies to the damped harmonic oscil- 
lator equation. If we consider only the radial motion (for a 2D 
disc), we have Z = 0, and 



dv,. v~ 



dr R So R^ 
R 



(11) 



dr 



Si 



Even for discs in two dimensions, Eq. (fTTT l is not analyti- 
cally tractable. However, as some of the parameters involved 
in the equation are small, approximations of the solution can 
be found by performing perturbative expansions. So me of the 
classical results detailed below have been studied in and 
lNSH86l but are here properly justified. The principle of those 
expansions is described on Fig.[T] At f = 0, r = ro which implies 
that R{T = 0) - 1. Because of gas drag, a grain dissipates both 
its energy and angular momentum and therefore, experiences a 
radial inward motion, i.e. R < I. The first parameter with respect 
to which a perturbative expansion can be performed is rjo (linked 
to the pressure gradient by Eq. ( IC.3b ) as it takes values of app rox- 
imately 10^^-10"^ in real protoplanetary discs (see iNSH86i) . and 
thus 770 «; 1 . We consider that this inequality also implies that 



(12) 



This inequality is always verified when <? < 1 and thus ap- 
plies to observed discs (see Sect. |6]l. For ^ > 1, there is a region 
where this inequality is not verified. However, in this case, the 
pressure gradient has the same order of magnitude as the gravity 
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Equations 
of motion 



Ts « Tl< 

If p=0: W77 small grains 



S„ RP » 1 



T|„ R-'*i « 1 



Tl< « Ts 

If p=0: W77 large grains 



1o 



R-fl*i « 1 



NSH 
86 



A-mode 



, R-fl*i « 1 



S„ RP » 1 



B-mode 



Fig. 1. Principle of the various perturbative expansions for the 
grain radial motion. Expanding first with respect to the small 
pressure gradient {rjoR^''^^ 1) leads to NSH86 equations. 
Expanding first with respect to the grain sizes {S^R'' <^ I oi 
SqR'' » 1) leads to W77 expressions for the particular case 
p = 0. Combining both leads to the A- and B- mode, respec- 
tively for small and large grains. 



of the central star and the model of a power-law profile for the ra- 
dial temperature is not accurate enough to model realistic discs. 
We thus conside r that for real discs Eq. ( fTSl i is always justified. 
Then, following lNSH86l we consider the system of equations 
given by Eq. (fTTT i. We set 



(13) 



and look at the orders (9(1), Oirjo),... of the expansion - see 
Appendix|E] We find that: 

iSoR"-^- (i - y/i-mR-"^') , 

Vr = 770 VH + O [ro] = I+RIPSI ^ ^ y^''> • 



The pressure gradient term has been retained to keep the gen- 
erality, however since we assume that 770 1, we equivalently 
have 

T 

R 



1 

R ' 



(15) 



Thus, to order O (/70), 



V,. = -- 



J]oSoR"- 



1 



-R^i'Sl 



or equivalently, using Eqs. (IC.4l i and (IC.ll l. 



Vr = 



(16) 



(17) 



SqR'' is the dimensionless expression of the ratio fs/fk- Eq- (HSJ 
shows that R^^Sl <sc 1 or R^i'S^ » 1, and thus RI'Sq <s; 1 
or R^So » 1, resulting in asymptotic behaviours for the radial 
grain motion. These asymptotic regimes were first described by 
|W77| for the particular case p - 0. They correspond physically 
to two limiting cases: where the gas drag dominates, which we 



call the A-mode, and where gravity dominates, which we call 
the B-mode. In the next sections, we study and describe these 
two so-called "regimes of migration" or "modes of migration" 
before treating the global evolution of grains given by Eq. ( fTSI l. 

3.1. A-mode (Radial differential migration) 

The A-mode corresponds to the regime R''So <sc 1 (or equiva- 
lently fs/fk ^ 1)- In the A-mode, Eq. ( fT6b reduces to 



dr dR „ „ „. 1 



or equivalently 



Vr = 



rcl d InP fs 
v'k dr fk ' 



(18) 



(19) 



where the 0{RPS{)) has been neglected. In this mode of migra- 
tion, the stopping time is much smaller than the Keplerian time 
scale. Considering one grain's orbit around the central star, its or- 
bital velocity is forced by the gas drag to become sub-Keplerian 
in just a few stopping times, i.e. almost instantaneously. Thus, 
the centrifugal acceleration is not efficient enough to counter- 
balance the gravitational attraction of the central star, and the 
grain feels an inward radial differential acceleration. The gas 
drag counterbalances this radial motion and the grain reaches a 
local limit velocity in a few stopping times. We call the phys- 
ical process of the A-mode of migration "Radial Differential 
Migration". 

The A-mode of migration originates first from a perturba- 
tive expansion for 770 <s; 1 (rigorously for 770^ ^ 1) and 
second from a perturbative expansion for So ^ \ (rigorously 
for SoRP <*: 1). Formally, we have performed: lim lim [...]. 

■?ri<< l ';o«i 

Historically, the A-mode had been derived by |W77| to explain 
the radial motion of small grains. In his study, he neglected the 
radial dependence of the stopping time and assumed Sq ^ I 
(this approximation also implies that R^''S^-^ <*: I, as R < 1, see 
Appendix|Fl). 

It is straightforward to integrate the differential equation 
Eq. ( fTSl ) by separating the R and T variables. Noting that 
R (T ^0) ^ l,we have: 



(14) - If-;' + ^+ 5 ^ 0: 



R -■ 
T 



+ 9 + 2 I 0^ 

1 - R-P+l+l 



(20) 



[-p + q+\^riQSo 



- \i-p + q + 



0: 



R = e"'""^°^ 
j_ In(^) 
770^0 ' 



(21) 



The outcome of the dust radial motion comes from a com- 
petition between two effects. As the grain reaches smaller radii, 

(1) gas drag increases, which slows down the radial drift, and 

(2) the differential acceleration due to the pressure gradient in- 
creases which enhances the migration efficiency. Point (1) is re- 
lated to iopt, which scales as the surface density profile, while the 
acceleration due to the pressure gradient in (2) is related to the 
temperature profile (see Eq. ( IC.ll l). Depending on which process 
is dominant, the grain's dynamics can lead to two regimes: 
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- If -p + q + < 0: then as 



T{R) 



[-p + q+ ^)t]oS( 



^0{r-p^''^^-), (22) 



the time it takes the grain to reach smaller and smaller radii 
increases drastically, according to the diverging power-law. 
Importantly, this behaviour constitutes our definition of the 
grain "pile-up". Mathematically speaking, accretion onto the 
central star occurs in an infinite time, i.e. 



lim R = Q. 



(23) 



- If -p + q + ^ > Q: the grain is accreted onto the central star 
in a finite migration time given by 



1 



T]oSo[-p + q+ 2) 
which increases as and 770 decrease, so that 



Um = 0. 



(24) 



(25) 



The presence or absence of a physical grain pile-up is therefore 
demonstrated considering the asymptotic behaviour of R{T) at 
large times. It is important to realise that the pile-up is a cumu- 
lative effect that can not arise from velocities only (which how- 
ever provide qualitative information on the grain's motion) but 
can only be found by integrating the equation of motion. This 
rigorously allows us to distinguish two different behaviours for 
the outcome of the grain's radial motion, and thus two classes of 
discs with respect to the A- mode. 

3.2. B-mode (Drift forced by a resistive torque) 

Returning to Eq. ( fTSI l, the B-mode corresponds to the other 
asymptotic regime, wheieR''So » 1 (or equivalently fs/f^ » 1). 
In this case, 

rv^f^-f^^-^^i (26) 

dT So 



or equivalently 



rCj d InP fk 



(27) 



In this mode of migration, the stopping time is much larger than 
the Keplerian time scale. Hence, the orbital velocity of a grain 
around the central star is almost the Keplerian velocity. However, 
because of the pressure gradient, the gas orbits around the central 
star at a sub-Keplerian velocity. Thus, the azimuthal differential 
velocity between the gas and the grain generates an azimuthal 
drag force whose torque extracts angular momentum from the 
grain. Given that the Keplerian angular momentum increases 
with radius (/ oc -y7), this torque results in the inward migra- 
tion of the grain. We call the physical process of this B-mode of 
migration "Drift Forced by a Resistive Torque". 

As for the A-mode, the B-mode of migration can also be 
derived first from an expansi on in (S pR'') and then from an 
expansion in rjQ. Historically, |W77| found an expression while 
only assuming that » 1 since he considered a flat density 
profile. To find the expression derived bv IW77l for large grains, 
we must assume that SqR'' » 1. It is crucial to see that this 
expression does not imply that » 1 when p > and /? — » 0. 

It is straightforward to integrate the differential equation 
Eq. ( |26] | by separating the R and T variables. Noting that 
R(T ^Q)^ l,we have: 



1 

So 1 



■ RP+1+1 



(28) 



m [p + q+ff 



- If p + q ■ 



R^e 



10 



'-ln{R). 



(29) 



As for the A-mode, two kinds of behaviours appear, depend- 
ing on the p and q exponents; 

- If p + q + J < 0: The grain migrates inwards, piles-up in the 
disc's inner regions and falls onto the star in an infinite time; 



lim R^Q. 



(30) 



However, the negative exponents required to be in this 
regime do not correspond to physical discs. Therefore, the 
grain dynamics in the B-mode in real discs belong to the 
second case; 

If p + q + ^ > 0; The grain is accreted onto the central star in 
a finite time 



(31) 



j]o [p + q+\) 

which increases as increases and 770 decreases and so that 
hm R^O. (32) 

As for the A-mode, considering the limit of R{T) at large 
times also proves the existence of two classes of discs with re- 
spect to the B-mode of migration. The radial motion of a grain 
in the B-mode of migration is also driven by a competition be- 
tween the increase of both the drag and the acceleration due to 
the pressure gradient. However, for real discs, p > and ^ > 0, 
and therefore, p + q + ^ > 0. Grains migrating in B-mode in such 
discs fall in a finite time onto the central star 

3.3. Radial evolution and asymptotic behaviour of single 
grains 

As we have seen, the grains behaviour is divided into two asymp- 
totic regimes, called the A-mode and the B-mode, which come 
from two different physical origins. However, the two criteria 
determining if the grains are accreted onto the central object in a 
finite or infinite time differ for the A- and the B-mode. It is thus 
crucial to determine in which mode a dust grain ends its motion 
to predict if the grain is ultimately accreted or not. Returning 
again to Eq. (fT6b . we have 



d^ _ -rjoSoR''-''^'- 
dT ~ TTr^psT 



(33) 



We can separate the R and T variables and integrate to obtain an 
expression for T{Ry. 



T = 



1 



[Ta{R) + Ts(R)], 



(34) 
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I 

* 1 > 

1 R 

IfTlo/ V^X andT„\ 

Fig. 2. Radial evolution of grains in the {R,So) plane showing 
that a grain in the Epstein drag regime ends its radial motion in 
the A-mode. The solid curves represent R '' for various values of 
p, they separate the A-mode (below) from the B-mode (above) 
regions. The horizontal dashed lines show trajectories of grains 
as they migrate inwards from R = I. The shaded area is a for- 
bidden zone. 
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Fig. 3. Radial motion R (T) of dust grains in the Epstein regime 
for rjo = 10"^. varies from 10"'' to 10^. Top: p - 0, q = ^, 
here -p + q+ j > and the grain is accreted onto the central star 
in a finite time. Bottom: p - ^, q = j, here -p + q + \ < and 
the grain piles up and is consequently accreted onto the central 
star in an infinite time. 



where 



and 



Ta(R) 



Tb (R) 



' 1 - R-P+9+\ 

-p + q+\ 
-\n{R) 



if - p + q + \ i^Q 



(35) 



if 



■p + q + h 



0. 



'5— — — 1- ifp + ^+s^^O 



p + q+2 

-Sl\n{R) \fp + q + 



(36) 



0. 



Eq. ( l34l l provides the asymptotic behaviour of the grains at large 
times. Interestingly, as /9 > for realistic discs, the contribution 
of the B-mode becomes negligible when R <sc 5'^'^''. Hence, 
grains initially migrating in the A-mode stay in the A-mode, but 
grains initially migrating in the B-mode end their radial mo- 
tion in the A-mode. This behaviour is summarized on Fig. [2] 
and d etailed in Appendix |G] This result was not predicted by 
IW77i as he neglected the radial dependence of the stopping time. 
Mathematically spea king, it comes from the fact that the pertur- 
bative expansion of IW77l has been performed with respect to 
powers of and not powers of SqRP . 

To illustrate the radial motion of dust grains in protoplane- 
tary discs, we numerically integrate the equations of motion for 
different values of the parameters 770, ^o, p, q- We set /70 = 10"^ 
to mimic a realistic disc and vary the order of magnitude of 
from 10"'' to 10^ for two sets of {p,q) va lues. First, we choose 
{p = Q, q - |); according to the INSH86I expansion, the grains 
are accreted by the central star in a time T^. Second, we set 
(p - J, q - |); the grains fall onto the central star in an in- 
finite time from the same approximation. This set of {p,q) val- 
ues is taken to mimic discs profiles that are commonly used and 
for which -p + q + j can take a positive or a negative value. 
Consequently, we interpret the radial motion R(T) of dust grains 
plotted in Fig.|3] 

- The top panel of Fig.[3]shows the results for (p - 0, q = j): 
Grains fall onto the central star, initially in the A-mode for 



the small grains and in the B-mode for the large ones. The 
radial-drift process is long for small and large grains but is 
optimal for grains with S = 5 m = 1 for which the accretion 
time is = 1.6/770, or = 160 with 770 = 10"^. 
- The bottom panel of Fig. [3] shows the results for (p = |, 
q = |): In this case, the radial density profile is steep enough 
to ensure that the grains are not accreted onto the central star 
To reach a given radius (for example Rf = 0.1), the optimal 
size is Sm,f - 2.9 = <9(1) (see Eq. ( IG.lb ). Hence, in this case 
grains efficiently reach the disc inner regions without ever 
being accreted onto the central star. The transition from the 
B-mode to the A-mode (for which R oc T^'^ in this case) for 
the large grains is visible in this plot. 



4. Radial motion in the Stokes regime 

Radial migration of large particles occurs in the Stokes drag 
regime, which depends on the dynamical viscosity fi of the gas. 
For hydrogen molecules: 



5/77 y/n 
64cr, 



m 



(37) 



where 77? = 277?h = 3.347446922 x 10"^' kg and cr^ = 
2.367 X 10"'^ m^ is the mol ecular cross section of the molecule 
dChapman & Cowlinglll97Cil) . The kinematic viscosity v is then 
defined by ju = pgV and the gas collisional mean free path is 
given by 

A,^f-^. (38) 

We now generalise the procedure outlined in Sects. |2]and|3] 
to the three Stokes regimes of Eq. Q. Using the dimensionless 
coordinates described above, we have 



1^0 = 



5m y/n 
64^^ 



-CsO, 



(39) 



(40) 
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B-mode I 




Fig. 4. Radial evolution of grains in the (R,So) plane showing 
that a grain in the linear Stokes drag regime ends its radial mo- 

tion in the B-mode. The soUd curves represent R * for various 
values of q, they separate the A-mode (below) from the B-mode 
(above) regions. The horizontal dashed lines show trajectories of 
grains as they migrate inwards from R - I. Shaded area: forbid- 
den. 

where Cso is given by Eq. ( IB.23b . Thus, the expression of the 
kinematic viscosity v is 



1^0 

Pgo 



(41) 
(42) 



First, if < 1, the drag force is linear in v - Vg and thus has 
the same structure as for the Epstein regime. Comparing the ex- 
pressions of C (R, 0) for the Epstein and the linear Stokes regime 
(see Appendix let, all the results found for the radial motion in 
the Epstein regime can therefore be directly transposed by set- 
ting q' - q and p' - In this case, the grain radial motion 
does not depend on p anymore and the N^SH86 expansion of the 
radial motion for small pressure gradients provides (see Eq. (l33T l) 



d^ _ -riaSlR 
dT ~ 

These crucial results follow: 



4-1 



(43) 



- In the A-mode {S q ^ R * ), grains experience a pile-up and 
migrate onto the central star in an infinite time if -p' + q' + 



5 < 0, i.e. if ^ < -4 (which never occurs in real discs) 
- In the B-mode (5o » ■> 



star in an infinite time if p' + q' + \ < 0, i.e. if ^ < ^ 



), grains migrate onto the central 
i < 0, i.e. if ^ < |. 



Thus, similar to the Epstein regime, we derived one criterion 
for each mode and need to determine in which mode the grain 
ends its motion. For observed discs, ^ - 3 < (see Sect. |6]l, and 

4 

as the particle migrates inward, R becomes smaller than 
and grains end their radial motion in the B-mode (see Fig. |4|. 
This result is fundamentally different to the one we obtained for 
the Epstein regime. Indeed, for grains migrating in the A-mode 
in the Stokes regime at low Reynolds numbers, the criterion ob- 
tained for a pile-up in the A-mode is never satisfied for real discs. 
However, after migrating inside a critical radius, grains switch 



to the B-mode, for which the pile-up can potentially occur, de- 
pending on the value of q. The corollary is that in discs having 
q < \, i.e. a shallow enough temperature profile, large grains in 
the Stokes regime at small Reynolds numbers remain in the disc. 
Such a criterion is applicable for real protoplanetary discs. 

Second, if R„ > 800, the drag force is quadratic in v - Vg. 
Assuming that the radial motion is decoupled from the verti- 
cal motion, we perform the NSH expansion at small pressure 
gradient (cf. Eq. (fT3Tl). We find that whatever the integer j, 
{rjoR^py {\ - Vg) — > at the limit — > 0. This means that both 

Vr and Vg — s/l/R are flat functions as their Taylor series expan- 
sion equals zero at each order Consequently, they can not be de- 
termined by perturbation analysis. This property comes from the 
quadratic dependency of the drag with respect to the differential 
velocity and thus is not related to the grain size. Consequently, 
in this drag regime, the drag force is extremely efficient and the 
corrections to the Keplerian motion are negligible at every order 
of the perturbative expansion. The particles are very well cou- 
pled to the gas and do not migrate significantly. 

Third, for the intermediate case, we could not manage to per- 
form the expansion at small pressure gradients. However, we ex- 
pect an intermediate behaviour between the two Stokes regime 
at small and large Reynolds numbers. Consequently, if Rg > 1, 
the migration motion becomes less efficient as the drag force is 
no longer linear with respect to the differential velocity between 
the gas and the dust particles. Thus, the main constraint for the 
radial-drift barrier due to the Stokes drag comes from the low 
Reynolds number regime for which the migration motion is the 
most efficient. 

Finally, confusion often arises when defining the "radial-drift 
barrier" as the difficulty a grain has of "overcoming s = 5opt" 
(i.e.) reaching the B-mode. Indeed, as we have shown, grains 
can survive their migration motion in the Epstein regime when 
they are in the A-mode whenever -/9 + ^ + 1/2 < 0, and grains 
can start their migration motion in the B-mode but be accreted 
in a finite time if -p + q + 1/2 > 0. This study also shows for 
the Stokes regime that a grain ends its migration motion in the 
B-mode. However, as demonstrated in this work, the ability of 
the grain to overcome the radial drift barrier is only linked to 
the value of q.\f q > 2/3, the grain will be accreted onto the 
central star in a finite time, even if it has s > iopt- Thus, we 
would argue that the definition of the radial-drift barrier has to 
remain the ability of the grain to be accreted onto the central star 
or depleted from the disc within its lifetime. 



5. Limitations of the model 

We have demonstrated that the time it takes for grains to reach 
smaller and smaller radii increases dramatically under certain 
conditions. Specifically, if -p+q+j < (resp. q < |), grains ex- 
perience a pile-up in the Epstein (resp. Stokes) regime. However, 
the model developed for the radial evolution of dust grains in this 
paper remains simple in that we neglect several important physi- 
cal processes: turbulence, grain growth, collective motion of dust 
grains, dust feedback on the gas surface density and temperature 
profiles. We now discuss how those processes can modify the 
criteria derived above. 

1. The local press ure maxima created by turbulence 

dCuzzi et al.1 1200 iL |2008|) and flie c ollective effects due 
to the dust drag onto the gas phase (lYoudin & Goodm"anl 
l2005i) are known to slow down the dust particles. However, 
the efficiency of these processes — such as the non linearity 
of the streaming instability in global disc models and the 
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life time of the pressure maxima — in real discs remains 
difficult to quantify. Omitting these phenomena constitutes 
therefore an upper limit for the grain migration efficiency, 
which will be slowed by these additional processes. 
In this study, we assume that changing the dust distribution 
does not change the thermal profile of the disc. We also ne- 
glect the viscous evolution of the disc, assuming that the vis- 
cous timescales are larger than the characteristic timescales 
of the initial dust evolution. It implies that we assume that 
p is constant during the whole grains evolution. We can ex- 
pect that for long term evolution, the surface density profile 
will flatten, leading to a smaller value of p. This makes the 
Epstein criterion harder to be met, while the Stokes criterion 
is not affected. 

We have shown that even if the velocity of the grain's in- 
ward motion depends on their sizes, their outcome only de- 
pends on the surface density and temperature profiles. Thus, 
if we now consider growing (or fragmenting) grains, we ex- 
pect that (0 the intensity of the inward motion depends on 
the growth efficiency (this point will be discussed in detail 
in a forthcoming paper), but that (;/) the grains outcome re- 
mains determined by our criteria, regardless of the growth 
regime. 
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Fig. 5. Transition between the Epstein and the Stokes regimes 
in a protoplanetary disc of Mdisc - IO^^Mq extending from 
Tin = 10"^ AU to Tout = 10"^ AU for several values of p and q. 
Grains with sizes below (resp. above) than the curve experience 
the Epstein (resp. Stokes) drag regime. 



Following this discussion, our simple model deals with pro- 
cesses that are optimized for the grains to be depleted on the cen- 
tral object. Consequently, our Epstein and Stokes criteria for the 
radial-drift barrier constitute the least favourable limit for grain 
survival. Thus, we are confident when claiming that the radial- 
drift barrier does not occur in some classes of discs. One may 
however expect that more discs are retaining their grains due to 
the complementary processes mentioned above. 

It should be noted that the criterion -p + q + 1/2 quantifies 
the outcome of the radial-drift motion of the grains, but not their 
kinematics (which depends on the grain size, the grain density, 
etc.). Thus, we provide predictions for which discs will retain 
the largest mass of solid particles, but do not predict for which 
discs the radial migration to the inner di sc regions is t he fas test. 
Full simulations like those developed in lBrauer et al. I (l2008h are 
required to make predictions of the dust kinetics, even more so 
when the grain size evolution is driven by a complex model 
of growth and fragmentation. However, this study suggests that 
while complex simulations are useful to study the details of the 
dust dynamics, they are not required to determine the grains out- 
come. 

As a conclusion of this section, we have mentioned that the 
physics treated by our model is not exhaustive. In real discs, the 
limits -p + q + j =0 and q - ^ may be softened by the effects of 
additional physical phenomena. However, these neglected pro- 
cesses (such as turbulence and grain growth) tend to decrease 
the efficiency of the dust radial motion. Our predictions of when 
the "radial-draft barrier" does not occur therefore remain valid. 
Our model represents a powerful indicator for predicting the dust 
behaviour in discs with given power law profiles: we expect that 
(i) discs satisfying -p + q + j < retain their small grain popu- 
lation and that (ii) discs satisfying q < ^ keep their large solids. 
On the contrary, discs for which -p + q + j > (resp. ^ > |) 
likely lose their small (resp. large) particles. 



6. Application to observed discs and planet 
formation 

6.1. Validity of ttie criteria in real protoplanetary discs 

We now study how the criteria we derived can be applied when 
considering the physical evolution of grains in observed discs, 
with finite inner radii and finite lifetimes. The analytic expres- 
sions of the previous sections have been derived using dimen- 
sionless quantities. We now provide the physical timescales of 
the radial dust motion estimating the parameters involved in real 
protoplanetary discs. We consider in this section a typical CTTS 
disc, of mass Mjisc = 10"^ M© around a 1 Mq star, extending 
from Tin = 10"^ AU to rout = 10"' AU. The disc inner edge is 
chosen to correspond to the dust sublimation radius for a 1 Mq 
star, whereas its outer boundary is representative of the largest 
observed discs. Its vertical extent is set by the choice of the tem- 
perature scale. We take T(l AU) = 1 50 K, a typical value ob- 
tained bv lAndrews & Wi lliams' (2001) in their disc observations. 

The transition from the Epstein to the Stokes regime occurs 
when /Ig = ^, or 

and is represented in the (r, s) plane in Fig. |5] for this typical 
disc for different values of the surface density and temperature 
power-law exponents p and q. The Stokes regime is seen to apply 
to large bodies in the disc inner regions and for large values of 
both p and q. 

Disc lifetimes are generally thought to be a few Myr 
dHaisch et alj 120011: [Carpenter et alj l2005h . and thus we take 
fdisc ~ 10'' yr. For a grain starting at a distance ro from a 1 Mq 
star, the dimensionless value T^sc is therefore 

^disc 

''''~t^~ rl'' ~ [r„(AU)]3/2- ^^^^ 
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position: 



(AU) 



Fig, 6. Values of 5o as a function of grain size s and initial po- 
sition ro for a disc of mass Mdisc = 0.01 Mq extending from 
Tin = 10^2 AU to rout = 10^ AU, with grains of intrinsic density 
Pd - 1000 kgm"^ and rjo - 10"^, for p - and q = j (top) and 
P - \ and q - J (bottom). The thick line shows the limit be- 
tween grains that are accreted onto the star (fin < fdisc) and those 
that survive in the disc (fin > fdisc), i-C- the survival limit, in the 
Epstein regime. 



The dimensionless value of the dust disc inner radius (n, 
0.01 AU) for a grain starting at ro is 



ro ro(AU)' 



(46) 



The link between dimensionless and real grain sizes is made 
through the optimal size for radial migration. Considering first 
the Epstein regime, its midplane value has a radial dependence 
given by 



■^opt 



(^0) 



S(r) 



(47) 



The dimensionless size of a particle of size s starting its migra- 
tion at position ro is therefore 



V27rpd 



Sopt.O 



^0 



p 



(48) 



Values of 5o are plotted in the (ro, s) plane in Fig. |6] for a disc 
of total mass Mdisc = 0.01 Mq extending from rin - 10^' AU to 
'"out = 10^ AU, with grains of intrinsic density pd = 1000 kgm"^ 
and 770 = 10"^, for both p - and p - ^■ 

The dimensionless time Ti„ for a grain to reach Ri^ is given 
by Eq. ( l34l i. In combination with Eq. ( |46] |. this gives an expres- 
sion of Tin as a function of and ro, whereas Eq. ( |45] ) gives 
an expression of Tdisc as a function of ro. Equating them yields 
a second order equation in 5o as a function of ro, which can be 
solved to determine under which conditions a grain reaches the 
disc inner edge at the end of its lifetime. Using Eq. ( |48] ) gives the 
corresponding relationship between the grain size and its initial 



p + q+^ 



^0 '^0 



3/2 



{-p + q+{){p + q+^) 



(49) 



which is plotted as a thick line in Fig. |6l It separates the (ro, s) 
plane into regions in which grains reach ^n and leave the disc 
before it dissipates (fin < fdisc) or survive in the disc throughout 
its lifetime (fin > /disc)- We call this curve the survival limit. 

For -p + q + J > 0, illustrated by the case (p - 0,q - j). 
Fig. |6] shows as expected that most of the grains are lost dur- 
ing the disc lifetime. However, small and large grains initially 
in the outer disc survive, therefore even with this profile, the 
disc retains a fraction of its grain population before it dissipates. 
Moreover, one may expect growing grains that reach 5 = 1 to 
be inevitably accreted onto the star (unless the growth process 
is fast enough for grains to outgrow the fast m igrating sizes be- 
fore they leave the disc, see lLaibe et alj|2008l) . Such discs may 
not form planets, but their remaining dust content may still make 
them observable, although they would likely be faint. 

For -p + q + \ < 0, illustrated by the case (p - j,q - |), 
even though all grains fall on the star in an infinite time, some of 
them reach the disc's inner edge before it dissipates. On the con- 
trary, for ro > 350 AU, grains of all sizes remain in the disc. This 
is also the case for all (sub)micron-sized grains, whatever their 
initial location, as well as most of the grains up to 0. 1 mm. These 
grains likely make the disc bright and easy to observe, since they 
are the grains contributing most to the disc emission at IR and 
submm wavelengths. A large reservoir of grains is available to 
participate in the planet formation process, however a firm con- 
clusion on their survival as their size evolves would require in- 
corporating a treatment of grain growth, as discussed in Sect. |5] 

It should be noted that the disc used in these examples repre- 
sents a lower limit, as it is low mass and very extended. A more 
massive disc with a smaller outer radius would have a larger 
surface density, and the corresponding survival limit in Fig. |6] 
would be shifted vertically towards larger sizes and more and 
more grains of larger sizes would survive. 

The left panels of Fig. Q show the influence of the surface 
density and temperature profiles on the location and shape of the 
survival limit curve in the (ro, s) plane in the Epstein regime. 
Increasing p from to 2 shifts the curve towards smaller radii 
and larger grain sizes, as well as slightly tilts it clockwise. The 
outer disc region in which grains of all sizes survive extends in- 
wards, as well as the surviving population of small grains as the 
curve's lower branch shifts upwards and becomes flatter On the 
contrary, the steepening of the curve's upper branch, confining 
the population of surviving large solids to the disc outer regions, 
is less dramatic. Increasing q from 5: to | also tilts the curve 
clockwards, but shifts it towards larger radii and smaller grain 
sizes. However, its effect is more limited than that of changing p. 
A disc with a steeper surface density profile and a shallower tem- 
perature profile is therefore more efficient at retaining a larger 
quantity of small grains and up to larger sizes. Indeed, large p 
and small q values are required to meet the -p + q + ^ < 
criterion introduced in Sect. |3] 

An equation very similar to Eq. (|49] l can be obtained for the 
linear Stokes regime by replacing p and ^ by p' = ^ and 
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Fig. 7. Survival limits of grains for different values of p and q. Fig. 8. Isocontours of the survival time (i.e. the time needed to 

Left: Epstein regime, right: Stokes regime. Grains to the left of reach the disc inner edge at - 0.01 AU) of grains of size s 

the curves (fin < fdisc) are accreted onto the star whereas those to and initial position ro for q - j and different values of p. 
the right (?!„ > fdisc) survive in the disc. 



q' = q (since the equation of motion has the same structure 
for both drag regimes, see Sect. |4]i, and using the expression of 
■Sopt.o for that regime, given in Table ICTI Here iopt.o <^ T'^: the 
weak temperature dependence results in very little change for a 
large range of temperatures below or above our adopted value of 
r(l AU) = 150 K. The right panels of Fig.|2]show the survival 
limit in the Stokes regime for the different values of q (note that 
it no longer depends on p). When q increases, the curve's lower 
branch slides towards larger radii, making the survival of parti- 
cles in the Stokes regime less and less favourable. 

Given the form of Eq. ( |49l ) and its different expressions for 
each drag regime, it is not possible to compute analytically the 
survival limit for a grain transitioning from the Epstein to the 
Stokes regime as it migrates inwards. However, large values of p 
and q, for which the Stokes region is the largest, are not observed 
in protoplanetary discs (see Sect. 16.2b . and in practical cases the 
Stokes regime only applies to a small area of the (ro, s) plane. 
Small grains, which are detected in disc observations at IR and 
sub-millimetre (submm) wavelengths, are mostly subject to the 
Epstein drag. We therefore focus on that regime in the following. 

Equation ( l49b can give quantitative information about the 
outcome of the grain population. Replacing fjisc by any time t 
gives the location in the (ro, s) plane of grains reaching the disc 
inner edge (at r - rin) in that time f, which is therefore the sur- 
vival time of those grains. Its isocontours are shown in Fig. [8] 
for different values of p and for q - j. Only one value of q 
is shown as the q dependence is moderate, as can be seen from 
the left panels of Fig. |7] The fate of particular dust grains can 
easily be obtained from these figures. For example, in the con- 
text of disc observations, 1 mm grains initially at 100 AU fall 
on the 10^ yr contour for p - 0. Their survival time decreases 
to a few 10"* yr for p ~ 0.8, and increases again to values larger 
than 10^ yr as p increases. At an initial position of a few hundred 



AU, 1 mm grains survive longer than 10*" yr for any p, therefore 
long enough to contribute to the disc emission over its entire life- 
time. As noted above, such grains have longer survival times for 
higher p and lower q values. As another example, in the context 
of planetesimal formation, the survival time of a 1 m particle ini- 
tially at 1 AU is ~ 10^ yr for p - 0, decreases to ~ 10^ yr for 
p ~ 1, and increases again to ~ 10^-10"* yrfor p -2. The ability 
of such particles to remain for long enough in the disc to grow 
to larger sizes therefore strongly depends on the surface density 
profile. As a general rule, the survival of pre-planetesimals in the 
inner disc is favoured by small values of p. 

Similarly, replacing now rin with any radius r in Eq. ( |49l ) 
gives the locus in the (ro, s) plane of grains reaching that radius 
r at any time t. Alternatively, one can plot isocontours of the 
initial position ro of grains in the (r, s) plane at various times, 
thus showing the radial evolution of grains with the same initial 
position but different sizes. This is shown in Fig.|9]for {p - 0, 
q - \) and Fig.[TO]for (p - q - |). These plots make it easy 
to compare the radial evolution of any particle to any physical 
timescale of interest in the disc. In particular, they show that the 
disc still contains particles at all radii at the end of its evolution 
(t - 10^ yr). No grains are found to the right of the ro = 10^ AU 
contour, since this was the initial outer disc radius. In the disc 
with {p - Q,q - |), no grains between ~ 0.06 and ~ 0.2 mm 
remain, and grains of other sizes still present were initially in 
the outer disc. Given that the grains of sizes which contribute to 
IR and submm emission have come from a small fraction of the 
initial disc, this disc is likely faint. In the disc with {p - ^,q - 
|), only grains with i ~ 0.1 mm are absent from the very outer 
regions, and the observable grains come from a larger portion 
of the disc, likely making the disc brighter than in the previous 
case. 

As a conclusion, the analytic criteria derived above apply 
even when taking into account the finite lifetime (or inner radius) 
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Fig. 9. Time evolution (in nine snapshots from t - 10"^ to 10^ 
yr) of isocontours of the initial position of grains in the (r, s) 
plane for the disc with p = and q - The label for each 
contour can be deduced from its abscissa in the upper left panel 
at t - 10"^ yr. 
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Fig. 10. Same as Fig.|9]for p - ^ and q - \- 



of the disc. For most CTTS discs, the dust is in the Epstein drag 
regime (except for some extreme values for the grains sizes and 
discs profiles). Therefore, the grain's radial outcome is given by 
the value of —p + q + j. However, the transition between discs 
for which the radial-drift barrier occurs or not consists more of 




Fig. 11. Histogramm of the q parameters obtained from 
[Andrews & WilUamsl (I2005L l2007h data of 63 observed discs. 
The distribution is roughly comprised between 0.4 and 0.8 and, 
centred around 0.55. Approximately 90 % of the discs satisfy 
q < 2/3. 



a continuum around the value -p + q + j - than in the sharp 
transition predicted by the analytic model. Therefore, the radial 
motion of the grains has to be studied on a case-by-case basis 



for discs close to the transition 
shown above in this section. 



-p + q + I - 0, using the figures 



6.2. Constraining pliysical systems 

We now turn to observed discs and check if they meet our 
Epstein and Stokes criteria to determine whether the radial- 
drift barrier is constraining for planet formation. To estimate 
the values of the p and q exponents f or real discs, we use 
the re s ults o f disc modeling obtained by [Andrews & Williamsl 
(l2005ll200l from data on 63 discs in p Ophiuchi, Taurus and 
Aurigae. Using sub-millimetre fluxes measured at several wave- 
lengths, they fit a range of disc parameters assu ming a geomet- 
rically thin irradiated disc with opacities from iBeckwith et akl 
(Il990h . a gas-to-dust ratio of 100, a disc radius of 100 AU and 
zero disc inclination. The temperature exponent q is well con- 
strained by the observational data set: the histogram of most 
probable q values is shown in Fig. [TT] However, p is not well 
constrained and is usually assumed to be |. Very flat profiles 
with p < ^ and very steep profiles with P > I seem to be 
excluded (iDutrev e t al."1996 tlWilner et al.l | 200dl: iKita mura et al 
2002[ 
2007; 



Testi et alln20 03: Isell a et al.1 120091: [Andrews & William 
Andrews et al.l 120 09). We represent the disc distribution 



I 



modeled by Andrews & Williams from observations in the (p, q) 
diagram of Fig. [T2l the histogram of Fi g . fTTlis represented by the 
gray-shaded area and spread over a range of p values, taking into 
account that extreme values of p are less probable. The dashed 
line {-p + q + 1/2 - Q) represents the border between migration 
in an infinite time and accretion onto the central star for the A- 
mode of migration in the Epstein regime, while the thick dotted 
line [q = 3/2) represents that same border for the B-mode of mi- 
gration in the Stokes regime at low Reynolds number. The two 
black circles indicate the discs used as examples in Sect. l3.3l We 
have split the disc distribution in four regions in the (p, q) plane: 
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Fig. 12. Location of the different outcomes of radial migration in 
the (p,q) plane. Dashed (resp. dotted) line: limit between accre- 
tion without or with grains pile-up resulting in a finite or infinite 
time in the A-mode of the Epstein regime (resp. B-mode of the 
Stokes regime at small Reynolds numbers). Shaded area: loca- 
tion of observed discs. Black dots: discs used as examples in 
Sect. [331 



1 . region I: -p + q + <0 and q < both small and large 
grains experience the pile-up effect. Those discs are poten- 
tially observable and may favour planet formation. 

2. region 2: -p + q + j < and q > j: only small grains 
experience the pile-up eff'ect: even though such discs retain 
their small grains, the population of pre-planetesimals in the 
disc inner regions may efficiently be accreted onto the central 
star (at least until they reach the high-/?g Stokes regime). 

3. region 3: -p + q + \ > and q < |: if the pre-planetesimals 
can form before the entire distribution of small grains has 
been accreted onto the central object, they will remain in the 
disc and may constitute planet embryos. 

4. region 4-: -p + q + j > and q > ^: both small and large 
grains are accreted onto the central star. 

The Epstein criterion indicates that for q values in the 
range constrained by observations, discs which keep their small 
grain population, and are therefore likely to be bright in the 
IR and submm, should have p values approximately in the 
[1; ll range. This is ind eed what is found in most disc sur- 
veys dRicci et al.ll2010a llbl). On the contrary, smaller p values 
should correspond to discs which lose most of their small grains, 
an d are therefore more d ifficult to detect. This is what is found 
by [Andrews et al.l ( |2010|) . who pushed th eir previous observa- 
tions of the Ophiuchus star forming region dAndrews et al.l2009h 
down to fainter discs, finding for this new sample a median p 
value of 0.9, lower than for brighter discs. The criterion we de- 
rive in this paper for small grains in the Epstein regime pro- 
vides therefore the correct behaviour for explaining the range 
of p values of observed discs. However, this result has to be con- 
sidered carefully for two reasons. Firstly, the p and the q expo- 
nents determined from the observations have to be considered 
with their respective errors. Given these uncertainties, one may 
not be able to distinguish between a strict negative or positive 
value for -p + q + ^. Second, the boundary between the differ- 
ent zones of the {p, q) plane consists more of a continuum rather 
than a strict limit due to the finite lifetime/inner radii of the discs. 



The outcome of the grains may thus not be predicted when the 
value of -p + q + ^ is close to zero. 

Now turning to the Stokes criterion for large solids, Figs.fTTI 
and [12] show that the vast majority of observed protoplanetary 
discs have shallow temperature profiles (q < |) and are thus 
able to retain their population of pre-planetesimals. These discs 
are therefore relevant places to find evidence of planet forma- 
tion, provided small grains can efficiently grow to form pre- 
planetesimals. For the remainder of the disc population, the out- 
come of pre-planetesimals will likely depend on their ability to 
reach the high Reynolds number Stokes regime. However, the 
case of a steep radial temperature profile can be encountered 
in at least one particular situation: circumplanet ary discs which 
typica lly have temperature profiles with <? = 1 ( lAvliffe & Batd 
l2009h . In this environment, we predict from our Stokes criterion 
that planetesimals will be accreted onto the planet. The timescale 
of the planet formation by the core-accretion process, which usu- 
ally corresponds to the time req uired to release the gravitational 
energy of the accreted bodies (lPollacketalJll996h . is thus in- 
creased as the drag from the gas onto the planetesimals releases 
an additional thermal contribution. 

7. Conclusion and perspectives 

In this study, we have generalised the radial grain motion studies 
of |W77, NSH86 and YS02 for both the Epstein and the Stokes 
regimes, taking into account the effects of both the surface den- 
sity and temperature profiles in the disc. As observations do not 
provide direct information about the three dimensional structure 
of discs, radial profiles of surface density and temperature are of- 
ten described by power laws: E(r) - ^[/^'' and T{r) = Tq'' ''^ 
where both p and q take positive values. The radial dust be- 
haviour in those discs is governed by the competition between 
gravity and gas drag. The final outcome of the radial motion 
is set by two counterbalancing effects. First, the temperature 
increases when the radius decreases. Consequently, the devia- 
tion from the Keplerian velocity increases, which accelerates 
the dust's radial inward motion. At the same time, the surface 
density also increases, which increases the gas drag efficiency 
and slows down the dust motion. The competition between these 
two effects fixes the ultimate mode of migration of the grain (A- 
mode, where the drag dominates or B-mode, where the gravity 
dominates) and thus the final outcome for the dust motion. In this 
work, we have shown that it can be represented by an analytical 
criterion which depends on the drag regime. For the Epstein drag 
regime (in which the ultimate radial motion is in the A-mode), 
if -p + q + J > 0, the dust particle is accreted onto the central 
star in a finite time, and if -p + q + \ < 0, the grain pile-up 
results in an infinite accretion time and small dust grains remain 
in the disc. We have shown that, as expected, these conclusions 
are somewhat mitigated when taking into account the finite disc 
lifetime and finite disc size. However, the outcomes still remain 
similar: the bulk of the small grain population is lost to the star 
in the first case, whereas in the second case the disc keeps most 
of its small grains. A similar criterion is found for the Stokes 
regime at low Reynolds number: if q < |, the accretion time is 
infinite and large pre-planetesimals remain in the disc and can 
constitute the primary material for planet formation. However, 
the Stokes radial motion differs from the Epstein regime as the 
ultimate radial motion occurs in the B-mode. 

The observational consequence is that discs with a large pop- 
ulation of small grains should be strong emitters in the infrared 
and sub-millimetre and should be easier to observe, and that 
those having lost most of their small dust should be fainter and 
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harder to detect. This is indeed what is found: a large fraction 
of the observed discs have large p values whereas fainter discs 
tend to have lower p values (Andrews et al. 2010), in agreement 
with this Epstein criterion. In addition, most of the observed 
discs have ^ < |, allowin g them to retain also their large pre- 
planetesimals. As noted bv lRicci et al.l (l2010albl) . explaining the 
data requires a mechanism halting or slowing down the radial 
migration of dust grains. We show here that local pressure max- 
ima need not be invoked, but rather that the combination of ade- 
quate surface density and temperature profiles is sufficient. The 
p and q exponents used to reach our conclusions are of course 
strongly dependent on the model used to fit the data. However, 
even varying the fitting models, a large majority of discs still sat- 
isfy both conditions -p + q+\ < and q <\ - Consequently, the 
radial-drift barrier (or the so-called metre-size barrier when con- 
sidering an MMSN disc) does not appear to constitute a problem 
for planet formation for the discs that we do observe. 

Our conclusions presented in this study assumed that 
the grain size remains constant during its motion. However , 
observations tell us that grairi s do g row ('Testi et aP 2003 
Wilner etaD 120031: lApai et al.l 120051: |Lommen_etalJ Ilo 
2009j). Grain growth is studied in various theoretical studies 
dSchmitt et all ll997HStepinski & ValageasI Il997t ISuttner et alJ 



19991: iTanaka et all j2005l: iDuUemond & Dominikl 



GaraudI 



20071: iBrauer et all 12008 



2005 



Klahr & Bodenheimer 2006: 
Laibe et al. 2008; Birnstiel et al. 2009i). In a forthcoming paper, 
we will generalise the formalism developed here to explain the 
radial and vertical behaviour of growing dust grains. 
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Appendix A: Notations 

The notations and conventions used throughout this paper are 
summarized in TablelA.ll 
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Appendix B: disc structure 

B.1. Hydrostatic equilibrium 

At stationary equilibrium - 0), gas velocities Vg^, Vgs, Vg, 
and the gas density pg obey mass conservation and the Euler 
equation; 



with: 



r 



gMcf {r,z')d.fir,z')dz', 



1 



finz) = - - 



1 



(B.12) 



V 



r dr dz 

PgVg.Vvg = -VP + pgg. 



The solution of Eq. dB.ll ) requires; 

Vgr - Vgz - 0, 



(B.l) 



(B.2) 



To simplify Eq. ( IB. 101 ). we first use the following identity; 
1 dP 2^^"(^^^g) 



Pg dr ^ dr 
which becomes with Eq. dB.l II ); 

^dln(p) 



which ensures mass conservation. Projecting the Euler equation 



1 dP 

Pg dr 



dr 



■ dr ■ 



(B.13) 



(B.14) 



on e,; 



1 dP 

Pg dz 



QMz 



Noting that / (r, z = 0) = and integrating /i by parts provides; 



(z2 + r2) 



.2^3/2 ■ 



(B.3) 



h {r,z)^gMc. 



Assuming that; 

P^c^ (r,z)pg, 
and dividing both sides of Eq. (IB.3l l by c^, we have; 

'91n(csPg) _ QMz 

" "(^2 + ^2)3/2^2- 

Integrating Eq. ( IB.5I ) between and z provides: 

QMz'Az' 



Jo 



gMfir,z')d-jcf ir,z')dz'. 



(B.4) and; 



1 dP 
Pg d?- 



^dln(p) 



( QM 



dr 



QMr 



(z2 + r2) 



3/2 



(B.15) 
(B.16) 



(B.5) 



+clgMfdrci' -cl-^^j^ gMfd.cfdz' 



Pgir.z) = 



P{r) ^ 
cl{r,z) 



{r^+z'^f'\l(r,z')_ (B.6) 



Then, Eq. ( IB.lQI i becomes: 

,,2 



This expression can be simplified by the following approxima- 
tions; 

- In the special vertically isothermal case, where the sound 
speed depends only on the radial coordinate, Eq. (IB.6I 1 sim- 
pUfies to; 



r ' dr 
Noting that ; 



' gMfd.cfdz'-clgMfdrcf. 

r^ dr Jo 



(B.17) 



Pg (r, z) = Pg (r) e 



'cHr) 



1 



1 



r VPT? 

,2 



I- f' fd.cfdz:^ r^d^,cfdz'+ rfd,A'cfdz', 
dr Jo Jo or Jo 

(B.18) 

and integrating the last term of the right hand side of Eq. ( IB. 181 ) 
by parts provides; 



(B.7) r fdrd.cfdz'=fd,c;^- f d-jfd,.cfdz'. 

Jo Jo 



(B.19) 



- Further, assuming a thin disc ((^) « 1), a Taylor series Therefore, Eq. ( IB.17l i reduces to 

' "I 
(B.8) ' 



expansion of Eq. (IB.7b leads to; 

Pg (r,z) =pg(r)e"^ 

rc^ {r) 



QM ,dlnP 2 r~lr,^^ -2 n -2\^ / 

+ c: — — + 0Mc; ydrfdfC^ - d-> fdrc^ ^ j dz , 

° (B.20) 



with: 



H{r)^ 



Vk (r) 



(B.9) 



dr JO 
which can be more elegantly written as; 

+ gMc; f [V/ X .eedz'. (B.21) 



which is the classical scale height for vertically isothermal 
thin discs. 



v'g0 QM ^d\nP 

— = — + c,s - 
r r^ 



dr 



B.2. Azimuthai velocity 

The radial component of the Euler equations is given by: 



1 dP 



QMr 



Pg dr (^2 + ^2) 



3/2' 



where pg is given by Eq. 



Thus; 



Thus, the expression of the azimuthai velocity of such a disc can 
be separated in three terms, called the Keplerian, the pressure 
gradient and the baroclinic terms respectively. This last term is 
neglected in most studies. For a three dimensional disc, this term 
rigorously cancels for - constant. In this case, the flow is 
inviscid and derives from a potential, and isobars and isodensity 
surfaces coincide; thus, there is no source of vorticity and the 
(B.IO) azimuthai velocity depends only on the radial coordinate. This 
terms also cancels out for flat discs in two dimensions. If the 
disc is vertically isothermal, Eq. ( IB.21I ) becomes; 



Pgir,z) 



cl{r,z) 



(B.ll) 



QM ,dlnP 



1 



1 



dMc,. (B.22) 
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B.3. Radial profiles of surface density and temperature Appendix C: Dimensionless quantities and 

In this section, we consider that the disc surface density and the equations Of motion 

temperature (and thus the sound speed) depend only on the radial To highlight the important physical parameters involved, we set 



coordinate and are given by the following power-law profiles: 



•2(r) = 2„(^f''=i;,r-^ 



(B.23) 



For vertically isothermal thin discs, the vertical density is there- 
fore given by Eq. ( IB. 8b with the scale height given by Eq. ( IB. 91 ), 
which can be expressed as; 



with: 



//(r) = //o - 



(B.24) 



(B.25) 



The expression of pg compatible with the vertical hydro- 
static equilibrium and providing the power-law profile set by 
Eq. ( IB.23b is written as: 



Pg =p;o'-"e"^''^w. 



Indeed: 

3 2 



\J — cx 



(B.26) 



(B.27) 



- ^(V ■ 



Hence, with - y/lrip' ,.H'„ and x - p 



- + - 



Pg = 



^0 



-(p-?+l)el^ 



(B.28) 



which gives the correct surface density profile when integrated 
with respect to z- With this expression of pg, P (r) is given by: 



P(r) = c2(r)pg (r,z = 0)-4. 



^0 



(B.29) 



which ensures that: 



dlnP 
dr 



(B.30) 



and, using Eq. (IB.22I) . we find: 



^'gfl= J— -{P + ^ + ^\<-r-'' - &Mq{ - - 



(B.31) 



V'kO 



and introduce dimensionless quantities given by the 



following expressions: 



r/ro 


= R, 




Tin 


= R 


-q 


Z/Io 


= R 


-p 


I'k/vko 


= R 


1 

2 


///Ho 


= R 






= z. 




Vgfl/V'ko 






^ \ 




. Pg/PgO 


= R 





(C.l) 



with: 



(C.2) 



The dimensionless parameter t/o gives the order of magnitude of 
the relative discrepancy between the Keplerian motion and the 
gas azimuthal velocity. We note that: 



rcj dlnP 
Vk dr 



VkO 



-r]oR 



-9+1/2 



(C.3) 



Then, we set fko = sj and define 



r t 






= T 


tkO 




Vr 


dR 






VkO 


= dr='' 


Vfl 


de 






VkO 


= ^df-^' 


z 






= z 








dZ 




df ■ 



(C.4) 



Writing the coefficient C (r, z) of the drag force of Eq. ^ as: 



C(r,z) = CoC(/?,Z), 
and using dimensionless coordinates, we have: 
Fn/md C(R,Z) 



(C.5) 



Vko/fkO 



We also introduce: 



and 



(I'ko'koCo) ' 



-|v-vg|''(v-vg). (C.6) 



■Sopt.o - (vjl^fkoCo)' , 



■Sopt.O 



so that Eq. (IC.6I 1 becomes 
Fo/wd 



Vko/t] 



C{R,Z)^^ . 
-r — V - Vc 



kO 



r'(v-vg). 



(C.7) 
(C.8) 

(C.9) 
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vides 



Table C.l. Expressions of the coefficients y. A, iopt.o and C {R, Z) 
for different drag regimes. 



Drag regime 



Epstein 



Stokes 
(R, < 1) 



Stokes 

(1 < i?g < 800) 

Stokes 
(800 < R^) 



■5opt,0 



C(R,Z) 



1.6 0.4 



2pd 



18rov;;'^Eo 



,0.6 
V ' kO 



1 .32^01:0 
8pd ^/2nHo 



_ 2 z- 

1° 5jg 5 2«3-, 



dr 
dr 



R 



So" R^ 
(veo - ^) 



(E.l) 



At this stage, we do not know the order of . We show in 



the lemma of Appendix iDl that 



Oirio). Applying this 



lemma, we see that at the order (3(1), taking v,o = and 
vgo = (which ensures that ^ - O (rjo)) is a relevant 

solution for the equations of motion (which corresponds to 
circular Keplerian motion). Thus, 



Physically, ^opt.o corresponds to the grain size at which the drag 
stopping time equals to the Keplerian time at tq. In Table IC.ll 
we give the expressions of y. A, iopt.o, C (R, Z) for the Epstein 
and the three Stokes drag regimes. The dimensionless equations 
of motion for a dust grain are then: 



— - — + {R,Z)\\ ■ 
AT R Si 

dVfl VgVr 
+ + 

dr R 



R 



dFZ 1 dZ 
dr2 5;, dr 



{R^ + 4>oZ^f^ 
C(R,Z)\y-y,\^ 

z 



{R^ + 0oZ2) 



3/2 



= 

= 0. 
(C.IO) 



Appendix D: Lemma for the different expansions 

Lemma: Let x be either r or and ; the order of the perturbative 
expansion. If: 

- VrO = 0, and 

- v.ri can be written as a function of R (Vj, = / (R)) with / - 
(9 (1) of the expansion in 770, 

then, — ^ is of order O (no). 
dr ' 

Proof: 



dv,, dv,, dR 

— = — = Vrf (R) 

dr dR dr ■' 



m^rlf {R)+0{t]1)^0(7]o). 

(D.l) 



Appendix E: Epstein regime: perturbation analysis 

- Order (9(1): At this order of expansion, 770^ is negligible 
compared to ^. Thus, substituting Eq. ( fT3l l into Eq. (fTTT i pro- 



V90 = a/-- 



(E.2) 



- Order O (770): Applying the lemma in this order of expansion 



and noting that 

dvflo _ _ 1 
dr ~ 2 

Eq. (fTTT i becomes 



R-'l^v,. 



Vh + 



(E.3) 



= 



So 



= -^R-'l^v,., - 



Si 



(E.4) 

Solving the linear system Eq. ( IE.4b for (y^i , vei) provides 
1 ISoR"-^- (1 - ^Jl-7^oR-1^') 

^ '^'''^ Z~ (E.5) 

1 R--{l- yJl-jjoR-"^') 



^0 



1 +R2pSi 



In addition to the expression of v,. given in Sect. [3j we also note 
that 



+ lo^'ei + 



R rr^- 



(E.6) 



1+R2PSI 



which provides with Eq. ( fTSl ) 



vb= a - - 



7]oR-^^'- 



(E.7) 



R 2[l +R2PSI) 

Appendix F: Link with Iw77l 's original derivation 

Following historic reasoning for small grains (see 

Sect. 13. Il l, we perform a perturbative expansion of the radial 
equation of motion with the variable. We verify that taking 
the limit at s mall 770 provides the expression found for the A- 
mode in the INSH86I expansion. (Formally, we will show that 
lim lim [...] - lim lim [...] ). Hence, we set: 



V,.o +SoVr\ +Slvr2 +0(^S^-^, 

veo + Sovei + Slvg2 + 0(^S^•^ , 



(F.l) 
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where we have used for convenience the same formalism as for 
the expansion in 770 - see Eq. ( fT3] ) (noting of course that v repre- 
sents different functions). An important point is that the lemma 
of Appendix|D]holds when substituting to tjq. Therefore, sub- 
stituting Eq. dF.ll ) into Eq. ( fTTT ) provides the equations of motion 
for different orders of O (So)- 



- Order (9 Eq. (fTTT i provides (9(1) expressions for the ve- 

VrO = 0, 

(F.2) 



locities: 



In this order of expansion, the azimuthal velocity corre- 
sponds to the sub-Keplerian velocity of the gas. There is no 
radial motion. 
- Order (9(1): 




(F.3) 



- Order (9 (5 0): 



Vrl - 0, 



1 

^2 



R ' 

-^+T]oqR- 



9-1 



y^ - TjoR-iT^oR"-"^^- 
Finally, we obtain expressions for v,. and vg: 

2SoRP-^- {i - yji-mR-"^'^ + o{sl) , 



(F.4) 



ve 



--r^^R-c,+slR"^^- 



R ' 



(F.5) 



1 -^_+rioqR-''-' 



We now c ompar e the INSH86I expansion at R^S <K 1 (A- 
mode) and the IW77l expansion at 770 <s 1. 

- INSH86I: From Eqs. (O and Eq. ( lEJl i: 



V,. = -■ 



l+R^PSl 



Vg 



1 ^70 R-"^'- 2\ 



2 1 +SlR^P 



- IW77I small grains: From Eq. iF.Sh 



Vg = 



--ijoR-^+SlRP^^lrjo 



R ' 



"2 /, \+<^[^o) 



- mR- 



2 -H 



R^p-''-h+o[sl) + o{nfi. 

(F.7) 



Clearly, Eq s. (IF.6I 1 and ( IF.7l i are identical, demonstrating that the 
theories of |W77| and ,NSH86i are consistent. We also not e that 
if the simplification of Eq. (fTsT i is not perfor med, the two IW77l 
expansions directly appear as the expansion of lNSH86l in O (So) 

Now, in the case of large grains, we perform a perturbative 
expansion of the radial equation of motion with respect to ^ 
while assuming that SqR'' s> 1, and verify that taking the limit 
at small 770 provides the expression found for the B-mode in the 
INSH86I expansion. 

Taking the same precautions as for the previous expansions, 
we write: 

1 J 1 

Vr = Vri) + —Vrl +0\ 
■JO 



1 

Vg = Vflo + — v'ei +0 
■JO 



<f2 



(F8) 



Following the same method as for the small grain sizes expan- 
sion, we obtain: 



- Order (9(1): 



VrO = 0, 



Vgo 



(F9) 



It this order of expansion, the azimuthal velocity of the grain 
is the standard Keplerian velocity. 
- OrderO(J^): 



Vrl = -2 

vgi = 0. 



TjoR-i 

R ' 



R 



(FIO) 



The expansion at higher order is more complicated and will not 
be used for further developments. At the order O^-p^, we have 
for \'r and vg: 



2 



1 



rioR-" 



O 



'J_ 



(Fll) 



We now compare the expressions pro vided b y the|NSH8^ expan- 
sion at R^pSq » 1 (B-mode) and the I W77l expansion at 770 1 



(F.6) for the radial velocity: 
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with 




Fig. G.l. The discrepancy (bottom panel) between the exact mo- 
tion (top panel) and its NSH86 approximation (central panel) is 
negligible. This is illustrated plotting the radial motion of dust 
grains for - 10"^, t/q - 10"^ and for p = 0, ^ = 4 (solid) 



and p - j,q = ^ (dashed). Top: Ri (T), middle: R2 (T), bottom: 
relative difference {R2 (T) - Ri (T)) /Ri (T). 



- INSH86I: 



l+R'-i'Sl 



(F.12) 



JO ^ 
- IW77l large grains: 



.^R-P-'i-h+O 
So 





' 1 ^ 








\^0) 



(F.13) 



Once again, we show that the lW77l and rNSH86l theories are con- 
sistent. 



Appendix G: Asymptotic radial beliaviour of a 
single grain 

Noting R\ (T) the position of a grain integrated directly from 
the equation of motion (Eg. ( fTTT )) and R2 (T) the position inte- 
grated from the INSH86I approximation (Eq. (34\). we highlight 
(Fig. lG.ll ) that the dis crepancy between the motion from the ex- 
act equations and its INSH86I approximation is negligible (the 
relative error is lower than 10"^ for all the considered sizes). It is 
therefore justified to use the analytical results derived in Sect. [3] 
to interpret the grain behaviour Thus, from Eq. (l34l l. we see that 
the time for a grain starting at R = 1 to reach some final radius 
Rf is minimized for an optimal grain size S m,f given by 



kf{R{) 



-ln{Rf) 



ifp + ^+ i^tO 



(G.2) 



if p -h ^ + 



0. 



As shown in Eq. ( |34] |. the outcome of the grain radial motion 
depends on the value of -p + q + j. 



If -p + q + J < 0: 



lim R^O. 



(G.3) 



For such disc profiles, all grains pile-up and fall onto the 
central star in an infinite time. Indeed, the surface density 
profile given by p > ^ + i is steep enough to counterbalance 
the increase of the acceleration due to the pressure gradi- 
ent. Therefore, grains fall onto the central star in an infinite 
time, whatever their initial size. Such an evolution happens 
because the grains always end migrating in the A-mode when 
they reach the disc's inner regions. A crucial consequence is 
that grains are not depleted on the central star and therefore 
stay in the disc where they can potentially form planet em- 
bryos. 

If -p + q + ^ > 0: 

]im R^O, (G.4) 



where 



rjoSo 



1 



-P + q+2 P + <^l+2' 



(G.5) 



In this case, grains fall onto the central star in a finite time. 
The surface density profile given hy p < q+\i?, now too flat 
to counterbalance the increasing acceleration due to pressure 
gradient. We note that: 

- For small sizes (5o 1), - OiSorjo). 

- For large sizes (So » 1), = <^(^)- 

- reaches a minimal value for a size S m given by 



p + q + 



-p + q+i 



Therefore 

(S m) 



m^{p + q+\}j{-p + q + i) 



(G.6) 



(G.7) 



S'ni is of order unity and corresponds to an optimal size 
of migration. Values of Sm and rioT^aiSm) in the {p,q) 
plane are shown in Fig. lG.2l When 5 ^ 5ni, both the A- 
and B-modes contribute in an optimal way to the grains 
radial motion. 

In this case, grains can be efficiently accreted by the central 
star if 5 - - (9(1). Thus, they can not contribute to 
the formation of pre-planetesimals. This process is called the 
radial-drift barrier for planet formation. 



im,f 



p + q + 
-p + q+j 



xktiRi), 



(G.l) 
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